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Resume 

Interactions between pairs of magnetic domain 
walls (DW) and pinning by radial constrictions 
were studied in cylindrical nanowires with surface 
roughness. It was found that a radial constriction 
creates a symmetric pinning potential well, with a 
change of slope when the DW is situated outside 
the notch. Surface deformation induces an asym¬ 
metry in the pinning potential as well as dyna¬ 
mical pinning. The depinning fields of the domain 
walls were found generally to decrease with increa¬ 
sing surface roughness. A DW pinned at a radial 
constriction creates a pinning potential well for a 
free DW in a parallel wire. We determined that 
trapped bound DW states appear above the de¬ 
pinning threshold and that the surface roughness 
facilitates the trapped bound DW states in paral¬ 
lel wires. 

Pacs : 75.60.Ch, 75.78.Cd, 75.78.Fg 

Nowadays, tailoring materials at the nanoscale 
to fabricate devices with precise functionality is 
an intensive area of research. A multitude of ma¬ 
gnetic nanostructures are proposed as storage or 
logic devices. This devices are mainly based on 
moving magnetic domains separated by domain 
walls (DW). The motion of DWs can be control¬ 
led by magnetic field or electric current [1, 2, 3]. 

Typically, device applications use ferromagnetic 
nanostrips (flat nanowires). To reliably control the 
movement of DWs in such structures, the DWs are 
usually pinned at notches [4] although other geo¬ 
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metrical imperfection as kinks [5], antinotches [6] 
or protuberances [7] were investigated. Enhanced 
DW propagation due to the edge roughness [8, 9] 
or internal disorder[10] was equally predicted. The 
pinning potential was found to depend on the DW 
type (transverse or vortex) [11]. For a transverse 
wall, the constrictions create a single potential 
well. 

Cylindrical nanowires present an alternative to 
the nanostrips. Transverse walls (TDW) in cylin¬ 
drical nanowires present some interesting proper¬ 
ties as the Walker limit is completely suppressed 
and theirs velocity and precession speed depend li¬ 
nearly on the applied current [12], They propagate 
without deforming theirs internal structure and 
are mainly studied for microwave oscillators based 
on DW. Very few studies investigated the pinning 
and interaction of DWs in such geometry [7, 13]. 
Arrays of cylindrical nanowires can be easily fabri¬ 
cated using electrodeposition. One way to create 
geometrical defects in such wires is by diameter 
modulation[14], therefore rotationally symmetric 
constrictions can be fabricated. 

In this paper, the interaction between TDWs 
and the pinning of a TDW by a rotationally sym¬ 
metric constriction (notch) in ideal and realis¬ 
tic (surface rough) cylindrical ferromagnetic na¬ 
nowires are studied. Our results establish that a 
radial notch produce a symmetric pinning poten¬ 
tial with a change of slope when the DW is outside 
the notch. When surface roughness is present, it 
is found that the pinning potential becomes asym¬ 
metric and that the depinning field generally de¬ 
creases. Furthermore, a pinned DW at a notch pro¬ 
duces a pinning potential for a free DW situated 
in a parallel wire (magnetostatic pinning). Trap¬ 
ped bound states of two DWs in parallel wires can 
form above the depinning threshold, due to large 
oscillations in the spatial position. Surface defor¬ 
mation increases the probability of the trapped 
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DW states. This can strongly influence the control 
of DWs in cylindrical nanowire based devices. The 
influence of the temperature is discussed using the 
escape rate theory and the stochastic ID model. 
Our results show the importance of understanding 
the pinning mechanism and of mutual interaction 
between DWs in cylindrical nanowires. 

To investigate the dynamics of DWs, 3D mi- 
cromagnetic simulations were computed using the 
Nrnag package[15] which determine the spatial dis¬ 
tribution of the magnetization dynamics [16] sol¬ 
ving the Landau-Lifschitz-Gilbert equation. Fig. 1 
shows the simulated structures. Cylinders with a 
length of 900nm and 60nm in diameter were used. 
The cylinders were discretized into a mesh with 
a cell size of 4nm, inferior to the exchange length 
(~5nm for Ni). We consistently checked on a num¬ 
ber of test cases that smaller cell size discretization 
(3nm) does not influence the results presented be¬ 
low. We use the material parameters of Nickel : 
saturation magnetization /xoM s =0.6T, exchange 
stiffness A = 1.05 x 10 _11 J/m, 7=188.5GHz/T (g 
factor of 2.15), damping parameter a=0.015. No 
magnetocrystalline anisotropy is considered and 
the temperature is set to T = OK. The stable struc¬ 
ture in our nanowires is a V-shaped TDW separa¬ 
ting domains pointing along the cylinder axis (x 
axis). The average position of the DW center (x) 
is extracted for each applied field along with the 
azimuthal angle (ip) of magnetization in the yz 
plane. 

The constrictions were created as a sharp modu¬ 
lation of the nanowires diameter, having a length 
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Figure 1 - (Color online) Cylindrical nanowires 
with a radial constriction used in simulations : per¬ 
fect cylinder (a) and cylinder with radial surface 
roughness with amplitude of 2nm (b) and 4nm (c). 
The equilibrium position of a pinned DW is shown 
in each case. 


of 20nm (same as wall width) and lOnm in height. 
The surface roughness (Fig. 1(b), (c)) was crea¬ 
ted from a random profile along the cylinder axis, 
with a predefined amplitude deformation a and 
variable correlation length A, full revolved around 
x-axis therefore preserving the rotational symme¬ 
try. Two radial deformation were used : a soft 
roughness with cr=2nm and A between 3nm and 
lOnm, and a hard roughness with <r=4nm and A 
between 3 and 15nm. 

Single notch, perfect wire. First, the pinning 
of a DW by a single regular notch situated in 
the center of a cylindrical nanowire was investi¬ 
gated. To determine the pinning potential, a dc 
magnetic field was applied and the distance tra¬ 
veled by the DW (at equilibrium) was compu¬ 
ted. The normalized potential energy (exchange 
+ demagnetizing) is represented in Fig. 2(a). The 
potential profile is symmetric with a change of 
slope around x=10nm, deviating from the simple 
harmonic potential of symmetric notches in flat 
nanowires[17, 28]. When most of the DW is in¬ 
side the radial notch (lower elongation x<10nm), 
the DW position varies linearly with applied field 
(see details in Ref. [28]) and the pinning poten¬ 
tial is harmonic with a spring constant k] = 
3.21 x 10 _4 N/m. Above x=10nm, when most of 
the DW is outside the notch, the variation of DW 
position is still linear but with a higher spring 
constant kh, = 1-23 X 10 _3 N/m as the DW starts 
to deform[18]. The maximal equilibrium distance 
from the center of the notch is 30nm before de- 
pinning at H^ ep =±300 Oe. While it travels to the 
equilibrium position the DW also precesses around 
the cylinder axis as shown in Fig. 2(b) and (c). The 
relaxation time can be as long as 600ns. As the dc 
field increases, the DW tends quicker to the equi¬ 
librium position and the oscillation damps faster. 

It was found that the pinning profile of the ra¬ 
dial notch does not depend on the chirality of the 
TDW. It is identical for head-to-head and tail- 
to-tail DW pinned at the notch. The profile is 
also identical if the external perturbation is a dc 
spin polarized current (polarization P=0.7, non- 
adiabatic parameter £=0.05), although the ma¬ 
gnetization in the DW precesses in a steady state 
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Figure 2 - (Color online) Potential pinning energy of a concentric notch situated in an ideal (a) and 
rough nanowire (d). The insets show the DW position for an applied magnetic field of 2000e. The 
lines in (a) represent fits with harmonic potentials, in (d) are guide to the eyes, (b) and (c) show the 
DW position and DW precession angle for two dc applied fields for the ideal case and (e) and (f) for 
three dc applied fields for the rough wire. 


(no damping of ^ due to the spin torque com¬ 
pensation). In this case, the depinning current is 
2.4xl0 10 A/m 2 . 

The domain wall width suffers a small modifica¬ 
tion when submitted to the external field (inset of 
Fig. 2(a)). The width contracts by 2 nm just be¬ 
fore depinning from 20nm (without field) to 18nm 
suggesting a small deformation of the wall when 
pushed by the external perturbation. 

Single notch, rough wire. To account for im¬ 
perfections in the wire, we studied the effect of ra¬ 
dial surface roughness on the pinning by a notch of 
a DW. The difference between the two imperfect 
wires (as in Fig. 1) and the perfect one is repre¬ 
sented in Fig. 2(d). We observe that the soft sur¬ 
face roughness introduce a low asymmetry in the 
pinning potential of the notch. Without external 
perturbation, the DW is pinned inside the notch 
at a local pinning site, which is not in the center 
of the notch. The elongation from the equilibrium 
position varies with the local complex potential 
present. Even the soft roughness already modi¬ 
fies the depinning field which is slightly asymme¬ 
tric (+2700e,-2800e) and smaller than for perfect 


wires. 

When hard roughness is present, the asymme¬ 
try in the pinning potential can increase highly 
and the local pinning sites outside the notch can 
trap the DW. In Fig. 2(d), a case is shown of 
hard roughness that presents a large asymmetry. 
The equilibrium position of the DW without ap¬ 
plied field is situated almost at the left edge of 
the notch. The depinning field is highly asymme¬ 
tric, and lower in one direction than the perfect 
case. The asymmetry in the depinning field is of 
70Oe (+2300e,-3000e) and the relaxation time 
varies (Fig. 2(e) and (f)). The oscillation of the 
DW position, at higher elongation, can take lon¬ 
ger time to damp as the local potential depends on 
the local geometry of the roughness. At -3000e, 
just before depinning, the DW is resonantly ex¬ 
cited (frequency 20MHz) by the applied field as 
shown in panels (e) and (f). The amplitude of the 
spatial oscillation is 5nm, while the DW magneti¬ 
zation precesses around x axis. The resonance fre¬ 
quency can be measured experimentally and the 
local shape of the potential can be determined [20]. 
When the DW depins, it jumps from a local pin¬ 
ning position to another while the magnetization 
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can change its direction of precession as the DW 
can move backward. 

Several cases of hard roughness were stu¬ 
died (the most representative are compared in 
Ref. [28]). In each of them an asymmetry in respect 
with the center of the notch was obtained. The 
asymmetry in the pinning potential was found to 
depend mostly on the roughness inside the notch 
which creates a local equilibrium position which is 
not in the notch center. The asymmetry is more 
pronounced when the zero field equilibrium posi¬ 
tion is close to the notch edge. The depinning field 
was found to decrease, compared to the ideal case, 
at least for one field direction. Sometimes, when 
the DW is pinned at an edge of the notch (without 
applied field) the depinning field increases slightly 
(lOOe) for one direction. 

Parallel ideal wires. Next, the interaction bet¬ 
ween DWs situated in two separated nanowires 
was investigated : a free tail-to-tail DW (TTDW) 
situated in a perfect wire and a head-to-head DW 
(HHDW) pinned at a notch (as above) situated 
in the center of a lower wire (Fig. 3(a) inset). 
The wires are separated by a distance y=60nm 
or y=80nm. In this configuration, the free DW is 
trapped in the magnetostatic potential created by 
the pinned DW. This type of magnetostatic pin¬ 
ning in neighboring nanowires was demonstrated 
experimentally in nanostrips [21, 22, 23]. 

The interaction is detailed in Fig. 3 : in (a) 
the interaction energy is represented function of 
the horizontal distance x between the two DW 
when an axial static magnetic field is applied. The 
potential well is symmetric and the interaction 
energy between the DWs can be fitted[28] with 
the analytical expression of the multipole expan¬ 
sion of Ref. [19]. The free DW in the upper wire 
has its magnetization pinned in the direction of 
the magnetization of the pinned DW (y direction) 
when the magnetic field is inferior to the depin¬ 
ning field (H dep = ±105Oe at y=60nm and ±750e 
at y=80nm), even if initially its position oscillates 
with high amplitude (90nm). Above the first de¬ 
pinning field, the spatial oscillations of the free 
DW (panel (b)) are complemented by full preces¬ 
sion of the DW magnetization in the yz plane and 


the DW is no longer trapped by the pinned DW. 
However, at +1250e (y=60nm) there is a parti¬ 
cular oscillation of the two DWs that finalizes in 
a bound state, even though the free DW does not 
have its magnetization pinned initially in the y 
direction. As can be observed in Fig. 3(c), the 
free DW magnetization precesses four times be¬ 
fore being captured by the lower DW and have 
its magnetization pinned in the y direction. The 
particular initial oscillation frequency of the free 
and pinned DW (see also panel (e)) makes possible 
such a state in cylindrical nanowires. Above or 
below 1250e, the oscillation frequency of the two 
DW does not lead to a bound state. The energy 
of the bound state is twice the energy of the two 
DWs just before depinning (see panel (a)). Any 
DWs state with energy lower than 7eV and higher 
than 3.5eV (on the fitted curve) can be a poten¬ 
tial bound state, if the initial oscillation (spatial 
and angular) of the two DWs leads to a parallel 
orientation of the DWs magnetization at the clo¬ 
sest approach. This DW bound state above the 
depinning threshold does not depend on the sepa¬ 
ration distance between the wires. If the wires are 
separated by 80nm, a bound state at 80Oe appears 
above the the first depinning threshold (±750e). 

Parallel rough wires. To determine the in¬ 
fluence of surface roughness on the DWs interac¬ 
tions, we considered wires with hard surface de¬ 
formation as before (er=4nm, A=3-15nm). As ex¬ 
pected, even at zero applied field, the DWs re¬ 
lax at local potential minima giving different re¬ 
laxation positions than the ideal wires. The pin¬ 
ning potential well is wider and symmetric around 
£DW2=15nm, in the case presented in Fig. 3(d), 
for a separating distance of 60nm. The free DW 
can travel longer distances during the initial tran¬ 
sient oscillation before being captured by the lo¬ 
wer DW, and its equilibrium position can be as 
much as twice the maximum equilibrium position 
in smooth wires. The free DW is considered trap¬ 
ped by the lower pinned DW, when it stabilizes 
at an equilibrium position having its azimuthal 
angle in the same direction (y) as the lower DW. 
The height of the pinning potential well dimi¬ 
nishes from 3.5eV to 2.5eV before the first depin- 
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Figure 3 - (Color online) (a) Interaction energy, E, between a free TTDW and a pinned HHDW 
situated in parallel perfect wires as a function of longitudinal displacement. The inset shows the 
equilibrium position of DWs. (b) Oscillation of the position of upper DW (x-dwi) for different applied 
magnetic fields above the depinning threshold. The rotation of the magnetization of the two DWs (in 
the yz plane) and the azimuthal angle are shown in (c) for H=+1250e corresponding to the bound 
state. The squares depict the initial position and the triangles the final position, (d) Comparison of 
E for smooth and rough wires separated by 60nm. The lines are guides to the eyes. The inset shows 
the equilibrium position for rough wires, (e) Oscillation of DWs position at and above the depinning 
field in rough wires. The circles correspond to the instant of first parallel orientation of magnetization 
in both DW. (f) Oscillation of the upper free DW position in the rough wire for different applied 
magnetic fields which lead to trapped states. 

ning threshold. The depinnig field is asymmetric, 

H(fe p =-690e and +830e, and corresponds to the 
first large oscillation that results in precession as 
in the smooth case. In panel (e), we detail the os¬ 
cillations of the pinned and free DW at and above 
the pinning threshold. The filled and empty dots 
represent the time stamps where for the first time 
the magnetization of both DWs are parallel. At 
840e, the free DW oscillates faster and comes clo¬ 
ser (by 20nm) to the pinned DW at their first ali¬ 
gnment than at 830e and therefore synchronize af¬ 
terward. We deduce that the potential well length 
is of 120nm, if the DWs are separated by longer 
distances when they align the trapped state can 
no longer form. 

The surface deformation produce a notable dif¬ 
ference : in smooth wires, only a bound state was 


obtained above the depinning field; in rough wires, 
we obtained a large number of bound states above 
the depinnig threshold (Fig. 3(f)). Before being 
captured, the free DW does several large ampli¬ 
tude oscillations and precessions, increasing the 
number of large transient oscillations every 50e. 
The roughness facilitates longer transient oscilla¬ 
tions and the loss of the initial trapping (azimu¬ 
thal angle in y direction) at lower applied fields 
and therefore the precession of the free DW. But 
the roughness also facilitates the trapping of the 
free DW by the pinned DW (panel (f)). A trapped 
state is obtained even after seven precessions of 
the free DW for H=-100Oe. This is probably due 
to the fact that the free DW, in the rough wire, 
jumps between local potential minima having a 
lower precession speed (due to the lower applied 
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fields) than in the smooth case. At ±105Oe, the 
free DW always depins as the equilibrium position 
would be to far to synchronize its magnetization 
with the pinned DW. 

The free DW jumps sometimes larger distances 
with little variation of the external applied field : 
from llnm (zero field) to 30nm (-50e) and from 
-42.5 to -79.5nm (90Oe to lOOOe), which indi¬ 
cates deeper local potential minima around these 
positions. The local deformation (local potential 
landscape) modifies highly the energy of the trap¬ 
ped state, the state obtained at lOOOe (x£>iy 2 =- 
79.5nm) is a high energy state (9.5eV) while the 
state obtained at -lOOOe (x£>^2=71.5nm) is a low 
energy state (3.2eV). 

The influence of the temperature on the DW 
can be estimated considering that DW depinning 
is usually a thermally activated process. The effect 
of the temperature, on the escape of a DW from a 
single pinning potential, can be modeled using an 
Arrhenius-Neel type law[23, 25, 26] : 

T = T 0 e~ AV/kBT (1) 

where AV is the barrier height, and To the at¬ 
tempt frequency at zero temperature. The value 
of To is estimated in the range 10 7 -10 12 Hz[23, 25]. 
Considering an experimental time of 1ms as in 
Ref. [23], the magnitude of the potential barrier 
a DW could surmount varies from 0.23 to 0.51eV. 
For the case of the radial notch in the perfect wire, 
this corresponds to a change in the depinning field 
of 70e, as the height of the barrier at OK is 4.75eV. 
The temperature will have a stronger impact on 
the depinning field of the rough wires, of 250e in 
one direction as the potential is asymmetric and 
the potential barrier is weaker. For DWs in paral¬ 
lel perfect wires, the temperature will change the 
depinning field by less than lOOe, while in rough 
wires it will change it by 20-25Oe. 

The depinning probability of a DW from a pin¬ 
ning potential can be also estimated with the sto¬ 
chastic ID model[27]. Using the same parameters 
as the micromagnetic simulation[28], we obtain 
that at T=0K, the depinning occurs at 2990e as 
an abrupt jump, the problem being deterministic. 
At 300K, a non null probability of depinning exists 
already at 2850e (of 3%) and the probability be¬ 


comes 100% at 2970e. The depinning probability 
varies on a range of lOOe, similar with the esti¬ 
mates above. As the potential barriers are high 
compared with the thermal energy, we do not ex¬ 
pect a dramatic change in the depinning held at 
300K. 

Cylindrical constrictions require large depin¬ 
ning held compared to the propagating helds of 
DWs. The depinning held of the captured free DW 
by a pinned DW situated in a parallel wire is of 
the same order as in nanostrips. However, as no 
Walker limit is present, trapping of the free DW 
can be realized above the depinning threshold at 
higher helds. The roughness facilitates the trap¬ 
ping, which can be used in design devices based 
on cylindrical nanowires. 

In conclusion, we have presented a systematic 
study of the interaction between TDWs in smooth 
and rough cylindrical nanowires. The static pro¬ 
perties of the pinning potential in a radial constric¬ 
tion have been analyzed. The potential well is 
symmetric and a change of slope is observed when 
the DW is situated outside the notch. The relaxa¬ 
tion times are longer than in nanostrips. We exa¬ 
mined the interaction between two DWs placed in 
parallel wires and observed that bound states can 
exists above an initial depinning held, due to the 
large oscillation of the DW positions. 

We also investigated the case of more realis¬ 
tic wires, where surface roughness is present. The 
deformation amplitude can change drastically the 
form of the potential well and the depinning held 
which becomes asymmetric. The interaction bet¬ 
ween DWs becomes more complex due to the 
asymmetry, dynamical pinning and temperature. 
The roughness increases the probability to obtain 
trapped bound states after initial large transient 
oscillations. 
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